From the conserved Kuramoto-Sivashinsky 
equation to a coalescing particles model 
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Abstract 



The conserved Kuramoto-Sivashinsky (CKS) equation, dtu = —d xx (u+u xx +v%.), has 
recently been derived in the context of crystal growth, and it is also strictly related 
to a similar equation appearing, e.g., in sand-ripple dynamics. We show that this 
equation can be mapped into the motion of a system of particles with attractive 
interactions, decaying as the inverse of their distance. Particles represent vanishing 
regions of diverging curvature, joined by arcs of a single parabola, and coalesce upon 
encounter. The coalescing particles model is easier to simulate than the original 
CKS equation. The growing interparticle distance £ represents coarsening of the 
system, and we are able to establish firmly the scaling £(t) ~ y/i. We obtain its 
probability distribution function, g(£), numerically, and study it analytically within 
the hypothesis of uncorrelated intervals, finding an overestimate at large distances. 
Finally, we introduce a method based on coalescence waves which might be useful 
to gain better analytical insights into the model. 

o 

OO - 

^^ . Key words: Nonlinear dynamics, Coarsening, Instabilities 

PACS: 02.50.Ey, 05.45.-a, 05.70.Ln, 81.10.Aj 



1 Introduction 



The study of growth processes of crystal surfaces [T|2|3] has turned out to be 
a source of a variety of nonlinear dynamics. A first, general distinction should 
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be made between a crystal growing along a high symmetry orientation (e.g., 
the face (100) of iron) and one growing along a vicinal orientation (e.g., the 
face (119) of copper). In the former case, growth proceeds [Hf5] layer-by-layer 
via nucleation, aggregation of diffusing adatoms and coalescence of islands; in 
the latter, the surface is made up of a train of steps [6] which advance through 
the capture of diffusing adatoms (step- flow growth). 

The interest in the nonlinear dynamics of a crystal surface mainly comes from 
the observation that growth is often unstable [7J. Step-flow growth plays a 
special role because it allows for rigorous treatments and the original two- 
dimensional character of the growth may reduce to effective one- dimensional 
equations: an equation for the density of steps when steps keep straight and an 
equation for the step profile, when steps move in phase. The two cases occur 
during step bunching [7] and step meandering [8f9] . respectively. 

As for the nonlinear dynamics resulting from the instabilities, they may vary 
from spatio-temporal chaos [TD] to the formation of stable structures [TT] , from 
coarsening processes [12] due to phase instabilities [13] to diverging amplitude 
structures |14|. Recently, T. Frisch and A. Verga have found [i~5|Tl~6] that in 
special limito the profile u(x, t) of the wandering steps satisfies the equation 



d t u = -d xx (u + u xx + u 



now known as the conserved Kuramoto-Sivashinsky (CKS) equation. The "con- 
served" label becomes clear upon comparison with the standard Kuramoto- 
Sivashinsky eq., d t u = —d xx {u + u xx ) + u 2 . 

An equation similar to ([T]), with an extra propagative term ^u xxxi arises in 
step bunching dynamics with vanishing desorption [T7] and in the completely 
different domain of sand-ripple dynamics [18J. The propagative term can be 
removed using the transformation u — > u+ (j/2)x, which, however, introduces 
7— dependent boundary conditions [17] . Numerics and heuristic/similarity ar- 
guments give a coarsening pattern whose typical length scale grows as £ ~ t n , 
with a coarsening exponent n — |, both in the presence [TTfTB] and in the 
absence [TJ] of the propagative term. 



The linear stability spectrum of the CKS eq. has the form u = q 2 — q 4 . In 
many equations having the same uj(q), coarsening occurs because the branch 
of steady states has a wavelength which is an increasing function of the am- 
plitude. These steady states are unstable with respect to phase fluctuations 
and the profile evolves in time, keeping close to the stationary branch |19j . 



1 In the case of vanishing desorption and weak asymmetry in the attachment ki- 
netics to the steps. 



For the above reason, our first step will be to discuss the periodic stationary- 
states of the CKS eq. (Section [2]), showing they look like sequences of arcs 
of a universal parabola, connected by regions of diverging curvature (asymp- 
totically, angular points). Direct simulations of the CKS eq. [Tb\ show that: 
(a) during the dynamical evolution, the interface profile can be thought of as 
a superposition of parabolas, and (b) the typical size £ of parabolas (in the x 
direction) grows with time as £(t) — \fi. 

In order to understand better the dynamics, we have simplified the problem, 
starting from the observation that angular points can be seen as effective 
particles interacting through the connecting arcs of the parabolas. The corre- 
spondence of a Partial Differential Equation (PDE) with a system of particles 
is known for the deterministic Kardar-Parisi- Zhang equation [20J, where arcs 
of parabolas are separated by cusps. In that case, the PDE is linearly stable, 
which translates to parabolas of decreasing curvature. For the Burgers equa- 
tion [21], cusps are replaced by shock waves and arcs of parabolas by linear 
pieces of decreasing slope. 

In principle, the particles move in two dimensions (the xu plane) and each 
particle interacts with all other particles. This full description, if possible, 
would be exact. However, we keep things simpler by limiting to the horizontal 
motion and to nearest neighbour interactions. Therefore, in Section [3] we show 
that the dynamics of particles is described by the equations &i = (xj+i— Xj) _1 — 
(xi — Xj_i) _1 , where Xi is the coordinate of the i— th particle and x i+ i > Xj. 
When two particles collide (xj = x i+ i), they coalesce and the total number of 
particles decreases by one, thus leading to a coarsening process. 

Subsequently, we simulate the particle model (Section 13.11) finding the coars- 
ening law I(t) and the size distribution of interparticle distances, g(£/£(t)). 
A very crude approximation of the Fokker-Planck equation for the particle 
system (Section I3.2J) gives the correct expression for I(t), but overestimates 
g(s) at large s. Finally, we introduce the method of coalescence waves (Section 
13. 3p . finding some preliminary numerical results, which might guide a future, 
more rigorous analytical study. 



2 Steady states 



The steady states of the CKS equation ([T]) satisfy the second-order nonlinear 
differential equation u + u xx + u x = a + bx . Since the constant b must vanish in 
order to get bounded solutions, while the constant a can be trivially absorbed 
into a uniform shift of u(x) and be set to zero, the problem reduces to solving 



the differential equation 

u xx = -u-u 2 x . (2) 

This equation also gives the steady states of a different PDE, studied by 
Mikishev and Sivashinsky [22J. Therefore, we limit ourselves to just a few 
results that play a major role in what is to follow. 

Interpreting x in Eq. (|2j) as time, the equation corresponds to a harmonic 
oscillator subject to an external force proportional to the velocity squared. 
Deriving (jSJ) with respect to x and putting u xx = F, one obtains 

*i = L_ (3 ) 

dF 1 + IF ' v ; 

Assuming the initial conditions w(0) = A and u x (0) = 0, the solution is 

F=-u-«* = -i + (i-A) e 2 ^). (4) 

We then have the trajectories in the (w, u x )-phase space: 



u ,. 



2 _ I _„ -(I- A\JiA-«) 



2 



u-{\-A)e 2 ^ . (5) 



If A > | the force F is strictly negative and the trajectory is not limited. 
On the other hand, if A = A + , with < A + < | , we get periodic, bounded 
trajectories, which oscillate between A + and — A_. For A + = |, we obtain 
the separatrix w 2 = | — w, which corresponds to the parabolic trajectory 
u(x) — \ — (x — x ) 2 /4. See Fig. [I] for more details. 

It is useful to determine the amplitude A_ in the negative x direction, as 
a function of A + . For A + — > 0, it is easily found that A_ — > A + , while 
in the important limit A + — > |, A_ diverges logarithmically according to 
A_e~ 2A - ~ i - A+, i.e., A_ ss -| ln(l - 2A+) . 

In proximity of u = — A-, bounded trajectories have a minimum with a cur- 
vature that diverges as A_ —> oo, as shown by the expansion u = —A_ + e in 
Eq. (jSJ), which gives u ~ —A_ + ^(5x) 2 + (^[^^(^a;) 4 ] . This approximation 
is shown as thin full lines in Fig. [T] (small upper inset). The quadratic and 
quartic terms are of the same order when \Sx\ ~ l/y/A_, which sets the size 
of the high-curvature region. In fact, the slope u x (5x = l/\/A~Z) ~ \J A_ joins 
the corresponding slope of the limiting parabola u = ~ — x 2 /4, when u ~ —A_. 

Finally, in the large lower inset of Fig. [T]we plot the wavelength A of the steady 
states as a function of their amplitude A_. The full line, A = Ay/~A~Z, gives the 




Fig. 1. Periodic steady states configurations, corresponding to A_ = 1 (dotted), 
A_ = 2 (short dashed), A_ = 4 (dot dashed), and A_ = 8 (long dashed). The thick 
full line corresponds to the limiting parabola u = -^ ~ x 2 /A. Small upper inset: the 
minima of the configurations and the approximations u = —A- + A_x 2 /2 (thin full 
lines). Large lower inset: wavelength A as a function of the amplitude A_. The full 
line is the asymptotic relation A = 4y9L~ . 

analytical approximation valid for large A_. It can be determined from the 



asymptotic parabola u 



x 2 /4, imposing w(A/2) = — A_. 



2.1 Steady states and dynamics 



In the Introduction we have argued that steady states are important because 
dynamics proceeds by evolving along the family of steady states of increasing 
wavelength Ao In the case of the CKS equation special attention should be 
paid to the constant a and to the conserved character of Eq. ([1]). We now show 
how the conservation law fixes the value of a, as a function of A. We shall then 
argue that a is the time- dependent vertical shifting of the surface profile. 



On the one hand, the conserved dynamics requires that the spatial average 
(u(x,t)) is time independent; on the other hand, steady states u(x) found 
in the previous section have a non vanishing and A-dependent average value. 
Therefore, the family Ud(x) of steady states which is relevant to the dynamics 



2 This is both an observation (see the following discussion on Fig. [2]) and a conse- 
quence of the stability of steady states with respect to amplitude fluctuations, see 
Ref. [151. 



is Ud{x) = u(x) + a(A), where the constant a satisfies the condition o(A) = 
— (u(x)). For large A, the average value of u(x) can be safely determined by 
approximating it with the arc of the (separatrix) parabola, so that 

(«(*)> = t / dx[^-^)+o(X 2 )^-^ + o(X 2 ). ( 6 ) 



Therefore, for large A we get Ud(x) = u(x) + A 2 /48. 



1 1 dx {\- 


X 2 

"T 


-A/2 





3 The particles model 



The following approach is founded on two observations, the first based on 
theory, the second on numerics: (i) When increasing the wavelength A of steady 
states, u(x) tends to a sequence of arcs of the universal parabola u(x) = 
a— (x — x) 2 /4, connected by regions of diverging positive curvature whose size 
5 ~ 1/A is vanishing small. See Section [1] and Fig. [TJ (ii) Dynamics deforms 
the above picture, but interface profiles can still be thought of as a sequence of 
points (x n , y n ), with x n+ \ > x n and the following properties: u"(x n ) — > oo with 
increasing time, and between any pair x n ,x n+ i of consecutive points u(x) can 
be approximated as an arc of the universal parabola, with a and x determined 
by the conditions u(x n ) = y n and u(x n+ i) = y n +i- See Fig. 3 of Ref. [15] and 
our Fig. [2J 

Therefore points (x n ,y n ) (henceforth particles, full dots in Fig. [2]) define un- 
ambigosuly the full interface profile and their dynamics should be derivable 
from the CKS eq. The coarsening process occurs because bigger parabolas eat 
neighboring smaller ones. When one parabola disappears, two particles merge 
into one: it is a coalescence process. Because of the conservation of the or- 
der parameter, d(u)/dt = 0, particles do not move independently and their 
effective interaction is expected to be fairly complicated and long-range. 

Here we limit ourselves to a simple model where a single bigger parabola 
(-Ui(x) = A — (x + h) 2 /4) eats a smaller one (w 2 = B — (x — /i) 2 /4), see Fig. [3] 
(full line). The full configuration is defined by the three parameters A, B and 
Xq, the point where the two parabolas meet. Two conditions must be fulfilled: 
continuity in xq implies B = A — hxo and conservation d(u)/dt = implies 
A + B — ^ — ^. = 0. In this picture, the particle in x = Xq has neighbouring 
particles at distances £- = 2(h + Xo) on the left and £ + = 2(h — xo) on the 
right. 

The next step is to use the CKS eq. to evaluate d(u 2 )/dt: it is enough to 




Fig. 2. Empty circles: The profile emerging from the dynamical evolution of the 
CKS eq. has been obtained digitalizing a late time profile in Fig. 3 of Ref. |15j . The 
full line is a sequence of arcs of the same parabola (x and u scales are arbitrary). 
Big full dots are located at the joins of the arcs and represent the "particles". 

multiply both sides of Eq. (fTl) by u and to integrate. We get 



1 d(u 2 ) 

2 dt 



{<) ~ (u 2 xx ) - ~((ul) x ) 



(7) 



where the last term on the right vanishes because of periodic boundary condi- 
tions. The average values (u 2 ), (u 2 ) and (u xx ) are evaluated using the two-arcs 
approximation depicted in Fig. [3] and are therefore functions of xq only. Once 
we replace (d/dt) by xo(d/dxo), we obtain a differential equation for the posi- 
tion xq of the particle, 



X 



2(» - (uU) 
d{u 2 ) 

dx 



The evaluation of (u 2 ) and (u x ) is straightforward, because the region of the 
angular point can be safely neglected: 



(u 2 ) 
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-^ + function (h) 
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{u x ) 12 + 4 ■ 



(9) 
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Fig. 3. Two arcs of the universal parabola, iti(x) = A — ^(x + /i) 2 and 
U2(x) = B — j{x — h) 2 , join in x = xq. This condition and the conservation of 
the average value {u(x)) fix two of the three parameters A,B,xq. The dashed line 
corresponds to xq = 0, the symmetric (and unstable) configuration. 



As for the average value {n% x ), we may write {u 2 xx ) = \ + (wiLKng = \ + 
(— f 2 -) , where (■ ■ -) ang means the average on the angular point and Au x = 
:)+ — (u x .)- is the abrupt change of slope occurring through the angular 

point, on a distance of order 5. Since Au x = u' 2 (xo) — u[(x ) = — (x — h)/2 + 

(x + h)/2 = h, we get {u 2 xx ) = \ + A ~ ^, the term | being negligible for 

large A (which means large h and small <5). 



2ft 



For a periodic configuration of wavelength A, we have 5 ~ 1/A. In the config- 
uration of Fig. [3] we do not have a single A, but the two quantities £±. If we 
assume 5 ~ (i^ 1 + £Zj, we get 5 = 5oh/(h 2 — Xq), with 5o being determined 
by the condition that the right-hand-side of Eq. ([7j) vanishes for xq = 0: in 
fact, the particle has zero speed in the symmetric configuration xq = 0. So, 
we get (w 2 ) — (w 2 ^) = ^2 + ^ H — 5^ : -^ ^f. It is worth stressing that a 

completely different assumption, 5 - 
indistinguishable from this. 



r 2 

6 ' 

-1(1-1 



gives a result which is almost 



We now have to determine the derivative appearing at the denominator of 
Eq. (|Sj): d(u 2 )/dxo = x (h 2 — Xq)/3, so we finally get x = x /(h 2 — Xq). 
Expressing x in terms of the interparticle distances £± = 2(h^fx ), we finally 
get £+£. = A(h 2 - x 2 ), (£- - £+) = 4x , and 



dx 
~dT 



1 



:m 



This is one of the main results of the paper: it means that the conserved 
Kuramoto-Sivashinsky equation can be translated into the motion of a sys- 
tem of particles with attractive interactions, decaying as the inverse of their 
distance, and undergoing a coalescence process when they collide. If the coor- 
dinate of the i— th particle is Xi, we can write 

dxi 11 . , 

~r = i r-7 n- 12 

dt (x i+1 - Xi) {Xi - Xi-i) 

Alternatively, if £i = (x i+ i — x,-) is the distance between particles i and (i + 1), 
we get 

« ± + ±-l (13) 



dt 



'i+l K -i—\ t-i 



Since interparticle force decays as the inverse of the distance, the coarsening 
law, i.e. the time dependence of the average distance between particles, £(t), 
can be easily inferred from scaling considerations, 

f~i - m~st. (M) 



This result is in agreement with numerical simulations [15] of the CKS equa- 
tion and with numerical results of the particle model, discussed in the fol- 
lowing section. It will also be corroborated by the Fokker-Planck approach, 
Section l3~2l It is worth noting that the result £ ~ yt is not related to diffusion 
since the particles motion is strictly deterministic; it arises from the l/£ decay 
of the interparticle force with distance £. 



3. 1 Simulation of the particle model 



Rather than solving Eqs. (1121) for particle positions, we have solved Eqs. (1131) 
for the distances. Fig. H] shows that the expected law £{t) ~ t 1 / 2 is satisfied 
for very different initial conditions: a random distribution (full diamonds) 
and a slightly perturbed uniform distribution (empty circles). The asymptotic 
value of £(t) is the same for the two distributions, showing that the prefactor 
c in £{t) = CQ\/t only depends on the initial density. The two sets (circles 
and diamonds) are distinct at small t because the initial random distribution 
favours more coalescences at short times than the uniform one. 

Fig. [5] shows the normalized distribution g(s) of interparticle distances as a 
function of s = £/£. Again, widely different initial configurations produce the 
same asymptotic distribution (empty circles and diamonds). As for the limiting 
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Fig. 4. The average distance between particles as a function of time, for the same 
initial density but different starting configurations: random distribution of particles 
(full diamonds) and a slightly perturbed uniform configuration (empty circles). The 
asymptotic law (full line) is £(t) = coy/i, with the same cq for the two distributions. 



behaviors of g(s), they are plotted in the upper inset for small s and in the 
lower inset for large s. At small s we clearly have a power law distribution, 
g(s) ~ s a , with exponent a ~ 1.3. At large s, data seem to suggest a gaussian 
tail. 

The qualitative features of the distribution g(s) can be understood by compar- 
ison with the one- species diffusion-limited coalescence process on the line |23j . 
In our case, particles separated by a smaller gap than average close the gap 
at increasingly larger speeds. As a result, the system evolves effectively as 
if nearest particle pairs react instantaneously, before the rest of the sys- 
tem evolves (numerical simulations do confirm this intuitive notion). In the 
diffusion-limited coalescence process, a quasi-static approximation shows that 
particle pairs react at a time proportional to the gap between them. Thus, 
both processes favor faster reactions between nearest particle pairs (in our 
case, more aggressively so). The distribution g(s) for diffusion-limited coales- 
cence is known exactly, yielding g(s) ~ s for small s, and a gaussian tail for 
s large. It is not surprising that we find a similar g(s). The faster-than-linear 
behavior, g(s) ~ s a , in our case, is consistent with the faster reactions be- 
tween nearest particle pairs, leading to a faster depletion of the probability 
distribution function near the origin. 
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Fig. 5. The asymptotic distribution of interparticle distances (empty simbols) for 
two different initial distributions (full symbols): random (diamonds) distribution 
and slightly perturbed (circles) uniform distribution. The system length is L = 10 6 
and there are 10 5 particles at the begininng. The asymptotic distribution is shown 
at time t = 10 5 , and is averaged over hundreds of runs. Upper inset: g(£/£) at small 
distances, on a log-log scale. Lower inset: g{£/£) vs. (£/£) 2 at large distances, on a 
lin-log scale. 

3.2 The Fokker-Planck equation 



In this Section we use the Fokker-Planck (FP) equation for the distances £{, 
under the approximation of uncorrelated intervals. The same approach has 
already been used for models of particles interacting with a force decaying 
exponentially, see Refs. [21125], and we refer the reader to those papers for 
more details. 



If f(£) = l/£ is the force between two particles at distance £, Eq. (TT3|) can be 
written as 



d£- 



(15) 



and the FP equation for the probability p({£},t) to find a given distribution 
{£} at time t is 







k 



d£ k 



(16) 
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We are mainly interested in the time dependence of the average distance, 
£{t), and in the probability distribution for the distances, g(£,t), which is 
expected to have the asymptotic scaling form g(£,t) = n(t)g (£/£(t)j, with 
n(t) = l/£(t). In the approximation of uncorrelated intervals, we get the 
equation d t g(£,t) = —2-^ (f(£) — f(£))g(£,t) , which does not include the 
coalescence process, because the form d t g = —d^J implies the conservation law 
dt£ = 0. Details on how to include coalescence are very similar to published 
papers [24.25] for different f{£), so here we merely state the results: 

£(t)=£ t 1 / 2 , (17) 

~g(s) = ~g - 8 ——. (18) 

s z - as + c 2 

The expression for £(t) agrees with numerical results (Fig. H|) and with scaling 
considerations, Eq. (fT4l) . On the other hand, the results for g(s) do not match 
our numerical findings, discussed in Section I3.1[ indicating the importance of 
correlations, neglected by this approach. The evolution equation (Tl3|) clearly 
shows that adjacent intervals are strongly anticorrelated, as gaps grow (or 
shrink) on expense of their surroundings; larger than typical intervals are 
surrounded by smaller than typical ones, and vice versa. Thus, larger than 
typical intervals grow far less than the independent interval approximation 
would allow, as the neighbors they engulf are typically smaller than average. 
In this way the approximation overestimates the frequency of long intervals 
(.9( s ) ~ V s ; instead of a gaussian tail). Similarly, being surrounded by large 
intervals, short intervals contract faster than if they were surrounded by typical 
intervals, explaining the overestimate of their frequency {g{s) ~ s instead of 
s a ) by the approximation. 

3.3 Coalescence waves 



This final section is dedicated to preliminary results on the study of coalescence 
waves, which might be useful for a deeper understanding of numerical results 
for the particles model. 

Consider a uniform distribution of N particles, perturbed at a single point: 
£i = £ for i 7^ N/2 and i^/2 = £ + A. As the system evolves, according 
to Eq. (fiH|) . the perturbation propagates on both sides of its origin, leav- 
ing behind a trail of coalescence events. These coalescence events take place 
at roughly equally spaced locations and time intervals, defining an apparent 
"front" , see Fig. [6j The speed v of the front (the slope of the straight line in 
the figure) is independent of the sign or magnitude of the perturbation, but 
seems to depend only on the density of background particles, v = Ci/£, with 
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Fig. 6. Temporal and spatial location of coalescence events following a single pertur- 
bation located in x = 5000. The two sets of data (empty circles and full diamonds) 
refer to perturbations of different intensity and sign. 
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Fig. 7. Full squares: Temporal and spatial location of coalescence events follow- 
ing two simultaneous perturbations in x = 4000 and x = 6000. These results are 
compared to the separate application of each perturbation (empty symbols). 

C\ ~ 4. The frequency, or inverse time between two consecutive coalescence 
events is roughly C2/£ 2 , with c 2 ~ 2. Therefore the coalescence front leaves 
a diluted system behind, with density n — > n' = (1 — c 2 /ci)n, or about one 
half the original density. Following these events, a second front of coalescence 
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events sweeps through, this time at roughly half the previous speed, due to 
the reduced density. This second front, clearly visible in Fig. [6j might be fol- 
lowed yet by others, but their quality deteriorates fast as the background of 
remaining particles distorts away from the original homogeneous spread. 

Next, we ponder how two propagating fronts, originating from two distant 
perturbations, interact. In Fig. [7] we plot simulations results for this scenario. 
One can see that coalescence fronts propagating in opposite directions anni- 
hilate. We have confirmed that annihilation takes place even when the two 
perturbations are started at different times. The rules for front propagation 
and interaction seem very simple, and give us hope that they might prove 
useful in shedding light on the kinetics of initially disordered particle systems. 



4 Summary 



The initial motivation of our work was to study the so called Conserved 
Kuramoto-Sivashinsky equation, Eq. (TO , whose stability linear spectrum is 
uj = q 2 — q 4 and whose numerical integration [T5] gives a coarsening process 
with an exponent n ~ \. The analysis of steady states (Sec. [2]) tells us that 
stationary configurations have the form Ud(x) = u(x) + a, where u(x) is a 
A-periodic function satisfying the differential equation u xx = — w 2 , — u and 
a = a(X) ~ — A 2 /48 is a constant fixed by the condition (ud(x)) = 0. 

Numerics and theoretical background suggest that the function u(x, t) evolving 
according to the CKS eq. keeps close to steady states. More precisely (Fig. [2]), 
u(x,t) appears to be similar to a continuous piecewise function, where each 
piece is a portion of the universal parabola y(x) = a — (x — x) 2 /4. Actu- 
ally, connecting points are vanishing regions of diverging positive curvature in 
the surface model. We have shown (Sec. [3]) that these angular points can be 
thought of as effective particles and we derived the equations governing their 
dynamics, Eqs. (I12lll3p . 



Simulating the particles model is definitely easier than simulating the interface 
model and we have obtained (Sec. 13.11) the coarsening law and the distribution 
of interparticle distances. We also expect that future analytical treatment 
would sooner be addressed to the particles model and we suggest two main 
directions: first, using the Fokker-Plank equation beyond the uncorrelated 
intervals approximation, used in Sec. 13. 2( second, using the coalescence waves 
method introduced in Sec. 13.31 
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